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Abstract 


Techniques were developed and described for performing three-dimensional finite element 
analysis of plain weave composites. This paper emphasizes aspects of the analysis which are 
different from analysis of traditional laminated composites, such as the mesh generation and 
representative unit cells. The analysis was used to study several different variations of plain 
weaves which illustrate the effects of tow waviness on composite moduli, Poisson’s ratios, and 
internal strain distributions. In-plane moduli decreased almost linearly with increasing tow 
waviness. The tow waviness was shown to cause large normal and shear strain concentrations 
in composites subjected to uniaxial load. These strain concentrations may lead to earlier 
damage initiation than occurs in traditional cross-ply laminates. 

Introduction 

Traditionally, advanced composite structures have been fabricated from tape prepreg, which 
was stacked to form a laminate. This type of construction tends to give optimal in-plane stiffness 
and strength. Since the primary loads usually are in-plane, this fabrication procedure appeared 
logical. However, there are at least two reasons why the usual laminated construction may not 
be best. First, secondary loads due to load path eccentricities, impact, or local buckling can 
sometimes dominate the failure initiation because of the low through-thickness strength of 
traditional laminates. Second, for thick laminates there are many laminae which must be 
assembled. This results in tedious labor with many opportunities for mistakes in orienting the 
laminae. 

Weaving is an alternate fabrication technique which has received considerable attention 
recently. The interlacing of fiber bundles in woven composites increases out-of-plane strength. 
Woven mats are thicker than a traditional lamina, hence fabrication of thick composites is less 
labor intensive and less prone to assembly error. These enhanced properties are obtained at 
the expense of some in-plane stiffness and strength. How much stiffness and strength is lost 
depends on the weave architecture. Because of the immense variety of possible weaves, it is not 
practical to determine optimal weave architecture through tests alone. Analytical models are 
needed whicli can predict the effect of various weave parameters on the mechanical properties. 

Most of the analytical models which currently exist for woven composites were developed 
for prediction of moduli (eg. ref. 1, 2, 3). These models are based on many simplifying 
assumptions, similar to those found in classical laminate theory, which may be appropriate 
for moduli prediction, but which preempt the extension of the model for strength prediction. 
Models developed for stress analysis have generally been quasi-three-dimensional (Q3D) (ref. 4, 
5). Since a Q3D analysis only models a single representative plane, little of the 3-D character 
of a woven composite is included. Ref. 2 includes some 3-D analysis, but the model was far 
too crude to permit stress analysis. Ref. 6 used a more refined 3-D finite element model than 
that in ref. 2, but the model was still probably not sufficiently refined for stress analysis (No 
stress distributions were reported.). There has been no detailed 3-D analysis of the stresses or 
strains in woven composites. 

This paper has two objectives. The first objective is to describe a refined finite element 
based 3-D analysis of woven composites. Figure 1 shows a schematic of the repeating unit 
for a single mat of a plain weave composite, which is the particular woven form considered in 
this paper. The grid which overlays the tows in fig. 1 is there to clarify the geometry of the 
tows. The discussion of the analysis will emphasize aspects of the analysis which would be new 
to persons who have only analyzed conventional laminated composites. The second objective 
is to present a few results which illustrate the effect of t ow waviness on the effective moduli, 
Poisson’s ratios, and the internal strain distributions. 
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Constitutive coefficients 

6x6 constitutive matrix 

Young’s moduli for orthotropic tow 

Average normalized Young’s moduli for woven composite 

Restraint forces in x-, y-, and z-directions, respectively, for 
load case i 

Shear moduli for orthotropic tow 
Average normalized shear moduli for woven composite 
Half-thickness of finite element model 
Half-lengths of wavy and straight parts of tow 
Applied load in the x-direction 
Displacements in x-, y-, and z-directions 
Specified displacements in x-, y-, and z-directions 
Half-width of finite element model ( W = l\ + I 2 ) 

Cartesian coordinates 
Normal strains 

Strains with respect to material coordinate system 

Poisson’s ratios for orthotropic tow 

Average normalized Poisson’s ratios for woven composite 


Analysis 

Three-dimensional finite element analysis was used to study the behavior of plain weave 
composites. The finite element method was selected because of its flexibility in modeling 
complex shapes, spatial variation of material properties, and arbitrary boundary conditions. 
For moduli calculations a crude model' can be used; for detailed stress analysis a much more 
refined model can be assembled, without any significant increase in professional manhours 
(compared to a crude model). In the following sections various aspects of the analysis will be 
discussed. First the configurations will be described, followed by sections on the finite element 
meshes and the material properties. 


Configurations 

The configurations analyzed consisted of mats of plain weave stacked to make a laminate. 
Figure 2 shows symmetric and unsymmetric stacking of the mats. The symmetry in fig. 2a is 
with respect to the interfaces between the mats. The geometry of an actual laminate is expected 
to be a complicated mix of these stackings. For this initial study the symmetric idealization in 
fig. 2a will be used. Such a laminate can be considered to be an assemblage of unit cells like 
that indicated in fig. 1. Because of symmetry within the unit cell, only one-fourth of the unit 
cell is actually modeled. Admittedly, this is a highly idealized plain weave composite. In actual 
composites the “synchrony” of the waviness in the mats could hardly be assured. Also, the 
interface between mats would not be planar. Furthermore, the shape of the tow cross-sections 
is likely to vary much more than is practical to model. In spite of the simplifications intrinsic 
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to the unit cell definition, this simplified cell definition is probably a reasonable starting point 
for detailed 3-B stress analysis. Fig. 3a shows a coarse finite element model of a plain weave. 
This model has a fairly small length of wavy tow (/ x ) compared to the length of straight tow 
(/ 2 ) Several different ratios of wavy to straight tows were considered. The ratio was varied by 
holding the length l\ constant and varying I 2 ■ The tow dimension II was also held constant. 
The term waviness ratio will be used as a measure of the fraction of tow which is inclined 
relative to the load direction. The waviness ratio is defined to be /j/(/ x + /j)- Note that tow 
waviness varies inversely with tow width for a plain weave; narrow tows result in very wavy 
weaves. 

Both extension and shear loadings were used. The boundary conditions for extension and 
shear were implemented quite differently, so they will be discussed separately. The dimensions 
and coordinate directions discussed below are defined in fig. 3a. 

The boundary conditions for extension loads involve constraints on both displacements and 
net normal force on the same plane. Because of symmetry, there are no shear forces on the 
planes. Imposing a specific displacement (zero or non-zero) is a simple matter in finite element 
analysis (eg. specifying u = 0 on x = —W). Imposing a constraint on both displacements am d 
forces on the same plane is much more complicated. For example, imposing the condition that 
v = constant on y = W and that the net restraint force is zero is not simple. Multi-point 
constraints could be imposed, but this is a complicated programming task and the number 
of nodes involved would result in a marked increase in the bandwidth of the equations. The 
approach taken herein is simple and does not affect the bandwidth. This approach, which is 
based on the principle of superposition, will be explained next by describing the steps required 
to obtain uniaxial loading of the unit cell in the x-direction. Because superposition is used, 
this approach is limited to linear analysis. 

The first step is to impose normal constraints on all six faces of the finite element model. 


constrain u on x — ±W 

constrain v on y = ±W (1) 

constrain w on z = ±H 

A displacement u = u Q is imposed on x = W and the other constrained displacements are 
set to zero. The normal constraint forces are calculated for the three planes x — W, y — W t 
and z = H and are defined to be F^ t Fy and F^ , where the superscript indicates that these are 
forces from loading case one. The subscript indicates the direction of the force (and implicitly 
the plane on which the force acts). Next, a displacement v =■ v 0 is imposed on y = W and the 
other constrained displacements are set to zero. The corresponding constraint forces Fy 
and F^ are determined. Finally, a displacement w — w 0 is imposed on z — H and the other 
constrained displacements are set to zero. The corresponding constraint forces F%, Fy and Fj* 
are determined. These nine constraint forces are used in the following equations. 

/'] + al'l + bf* = P 

Cj + aF% + bFy = 0 (2) 

F) + aFl + bFl = 0 

The unknowns are the load in the x-direction ,P, and the scaling coefficients a and 6. The 
last two equations express the condition that the net normal force on the y —W and the z H 
planes must be zero. Solving the last two equations yields 
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These values of a and b can be used in equation (2) to determine the load. The average 
normal strains and two of the average Poisson’s ratios are 
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2 W 
avo 

2 W 
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-h 

Sx 

e z 

Cx 


( 4 ) 


The effective Young’s modulus in the x-direction E x is calculated using the following energy 
balance equation 


\ Puo = \tixelVol (5) 

For the models considered in this paper one could have calculated E x based on just the 
the average stress on the plane x = W and the specified strain. In the more general case, the 
geometry might be such that there is no simple cross-sectional area which could be used to 
determine average stress. Such would be the case for the configuration in fig. 1 if the neat resin 
regions were left out. 

The remaining Poisson’s ratios and Young’s moduli can be determined in a similar manner. 

The exact boundary conditions for shear loading of an infinite array of unit cells require 
many multi-point constraints to impose anti-symmetry. Such boundary conditions permit the 
faces of the model to warp. The multi-point constraints complicate the analysis and significantly 
increase the computational cost. Hence, the approximate boundary conditions given below were 
used. 

for G xy * 

u — ay and v = ax on x = ±W and y = ±W (6a) 

for G yz : 

v — az and w ~ ay on y — ±W and z — db H (66) 

for G xz : 

u = az and w = ax on x = ±W and z = ±H (6c) 

These boundary conditions are very simple to impose. Consider the case of shear in the 
xy plane. Known non-zero displacements are prescribed on four faces. On z = ±H the 
displacement w is set to zero. This is posible because on the average, shear in the xy plane 
does not cause normal strain in the z-direction for the plain weave configuration studied herein. 
This is not generally true for textiles. For other material architectures a superposition procedure 
(like that described earlier) might be necessary. 



Finite Element Meshes 

Fig. 3 shows typical coarse and refined meshes for the plain weave composite. The coarse 
mesh had 595 nodes and 96 elements. The fine mesh had 3793 nodes and 768 elements. The 
elements were 20-node isoparametric hexahedrons. The meshes exhibit cyclic symmetry. In 
fact, a mesh is synthesized from the basic unit in fig. 4 rotated at 0°, 90°, 180°, and 270° about 
the z-axis. 

In order to minimize the number of elements required, compatibility of displacements was 
not completely maintained at the center of the model. Fig. 5 illustrates the incompatability. 
Elements A and B are the same as those so labeled in fig. 3a. They are shown separated 
in fig. 5 to show the interface between the two elements and to show element C. There 
is some incompatability in the displacements because node 9 is not connected to any node 
in element A. However, since the resin material (eg. element A) is much more compliant 
thau the tows, this approximation in the modeling is probably not significant. Even with the 
various simplifications, the internal shape of the mesh is fairly complicated. Work is needed to 
determine what refinements to the modeling would yield the greatest improvements in accuracy. 

Some of the meshes had a larger ratio of wavy tow to straight tow than other meshes. Fig. 6 
shows a mesh with a waviness ratio of .5, compared to a waviness ratio of .167 for the mesh in 
fig- 3. 

Material Properties 

The material properties chosen were for a hypothetical linear elastic graph ite/epoxy com- 
posite. The tow properties were selected to be the same as for unidirectional tape prepreg 
material. The resin was assumed to be isotropic. The table below summarizes the assumed 
material properties. 



Generic Graphite/ Epoxy* 

Neat Resin** 

27n 

13.42710 Pa 

.3452? 10 Pa 

E 22 

1.022710 Pa 

.3452? 10 Pa 

E 33 

1.022710 Pa 

.3452?10 Pa 

"12 

.3 

.35 

"23 

.49 

.35 

"13 

.3 

.35 

<7l2 

.5522710 Pa 

.1282710 Pa 

G23 

.3432710 Pa 

.1282710 Pa 

Gl3 

552JE710 Pa 

.1282710 Pa 


*ref. 7 
**ref. 8 


There are two angles which describe the orientation of a tow at any point. There is a 
primary rotation about the z-axis of either 0° or 90°. This angle is constant for a particular 
tow. There is also a secondary rotation about the x- or y-axis due to the waviness. The amount 
of rotation about the x- and y-axes varies spatially. This second rotation was determined based 
on the direction of the normal to the surface defined by four of the midside nodes, which are 
labeled a, b, c, and d in fig. 7. The surface normal direction was calculated as the cross product 
N = ca x db. Appendix A describes the material property transformations required because 
of these rotations. A right-handed coordinate system is assumed. Positive rotations are defined 
to be counterclockwise from the material axis to the global axis (or clockwise from the global 
axis to the material axis). 
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Results and Discussion 

Several configurations were analyzed to determine the effect of tow waviness on effective 
moduli, Poisson’s ratios, and internal strain distributions. The results of these analyses will be 
discussed in this section. 

Figs. 8a and 8b show the effect of tow waviness on the effective moduli and three of the 
Poisson’s ratios. The properties are normalized by those for a conventional (0/90) s laminate 
fabricated from tape material. For reference, the properties for a (0/90) s laminate are listed 
below. 


E x = 7.25£10Pa E z = 1.27E10 Pa 
G xy = .552£ r 10 Pa G yz = A26E10 Pa 
is xy “ .0424 ^ zy ~ .0^95 is x z ~~ -4o5 

Most of the results were generated using coarse models. A few results were generated using 
refined models. The differences in the moduli for the coarse and fine models were insignifcant, 
so only the results for the coarse mesh are shown in fig. 8a. There was more difference between 
the coarse and fine models for the Poisson’s ratios. Hence, results for both model refinements 
are shown in fig, 8b. The solid lines are the results for the coarse model. The solid circles are 
the results for the fine model. 

Fig. 8a shows that E x , E z and G xy decrease almost linearly with increased waviness. The 
waviness has the largest effect on E x . Some of the decrease can be attributed to the increase 
in resin content as waviness increases. The volume of neat resin pockets increases from 3.2 
percent for the case with a waviness ratio of .167 to 12.5 percent for the case with a waviness 
ratio of .5 . The transverse shear modulus G yz increases with increased waviness. Fig. 8b 
shows that the in-plane Poisson’s ratio v xy decreases with increased waviness. The other two 
Poisson’s ratios u xz and u zy increase with waviness. The variation is essentially linear. There 
is more difference between the results for the coarse and refined models for large waviness than 
for small waviness. Numerical values of the coarse model results in fig. 8 are tabulated in 
Table 1. 

No experiments were performed as part of this study, but some qualitative comparisons can 
be made with the results in ref. 9. Ref, 9 presents analytical and experimental in-plane modulus 
results for a plain weave, oxford weave, 5 harness satin weave, and an 8 harness satin weave. No 
results were presented for a (0/90/90/0) tape laminate, but the 8 harness satin results should 
have close to the same moduli, since it has low waviness. The loss in axial stiffness E x shown in 
fig. 8 agrees qualitatively with ref. 9. Fig. 8 shows a decrease in v xy with increased waviness, 
but the analytical results in ref. 9 predicts no change and the experiments showed an increase. 
Fig. 8 shows a decrease in G xy with increased waviness, which agrees with the experimental 
results in ref. 9, but which disagrees with the analytical results. There is obviously a need for 
further analytical and experimental work to clarify the source of these inconsistencies. 

Strain distributions were calculated for weaves with waviness ratios of .167 and .5. Refined 
finite element meshes were used. The procedure for smoothing the strains to obtain nodal 
strains is described in Appendix B. Because of limited space, only a few of the results will be 
presented here. Also, only uniaxial loading was considered for the results herein. The results 
will be presented in terms of isostrain contours. The strains are normalized by the magnitude of 
the applied axial strain. These strains are calculated relative to the material coordinate system, 
in which the xj-axis is along the tow direction, the X 2 -axis is in-plane and perpendicular to the 
xi-axis, and the X 3 -axis is perpendicular to the other two axes. Since no convergence study 
was performed, the results should be considered qualitative. 

For the case with a waviness ratio of .167, fig. 9 shows the 0° tow for which strain contours 
will be presented. Fig. 10 identifies the analogous area for the case with a waviness ratio of .5. 
The area in fig. 10 labeled ABCD is identical in geometry to the area labeled ABCD in fig. 9. 
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Hence, the only difference between the two models is the straight-tow region in the mesh with 
the smaller waviness ratio. 

Fig. 11 shows contours of constant t\ for the case with a wayiness ratio of .167. The largest 
strain concentration is at (x,y,z)=(0,0,0). There are strong, fairly complicated strain gradients 
near (x,y,z)= (0,0,0), but a few tow thicknesses away from the y=0 plane, the variation in strain 
on any y = y 0 plane is independent of y 0 . A two-dimensional analysis would probably do a 
reasonable job in the region away from the y = 0 plane, but obviously could not predict the 
peak strains, which occur near the origin. 

Figures 12a-12f show £j, s 2 , £3. £12. £23. and £ X3 contours for the region labeled ABCD 
in fig. 11. As just mentioned, the largest £ X occurs at the origin. The magnitude there is 
approximately 2.7 times the average strain for the unit cell. This strain concentration would be 
expected to cause fiber breakage at a lower global strain than would occur in a tape laminate. 
Much of the strain concentration is likely due to tapering of the tow in both the x- and y- 
directions to a zero thickness at the origin (see fig. 4). Herein lies a particularly sticky modeling 
problem. Even though the tow cross section may actually vary as indicated by the finite element 
model, the tapering is due to migration, not termination, of the individual fibers. That is, fibers 
near the y = 0 plane which are aligned parallel to the x-axis away from the origin are squeezed 
further away from the y — 0 plane as they approach the origin, resulting in a zero tow thickness 
at the origin. This fiber migration causes a very complicated variation in the local constitutive 
properties of the tow. There was no attempt in this paper to account for this. The straightening 
of the tow under tensile loads causes an increase in £ X at point E and a decrease in £ X at point F. 

Fig. 12b shows that the e 2 strains are small. The peak magnitude occurs on the y = 0 
plane, probably because that is the boundary of the tow, where there is an abrupt change in 
local properties. 

Fig. 12c shows that both the magnitude and variation of £3 are large. The presence of 
large £3 strains has been previously reported in ref. (4). In that study a quasi-3D analysis was 
used. Fig. 12c shows that the strain variation is nearly 2D in character, except where the peak 
occurs. The large magnitude of £3 might cause delamination. 

A uniaxially-loaded conventional (0/90) s laminate fabricated from tape prepreg has no shear 
strains away from free edges. Figs. 12d-12f show that a plain weave fabric has significant shear 
strains. The peak e l2 is about the same magnitude as the average axial strain. Particularly 
under fatigue loading, the £ 12 could lead to intra-tow cracking. The peak e 2 3 (fig. 12e) is also 
nearly as large as the average axial strain. The peak occurs along the tow boundary and might 
lead to some inter-tow cracking, particularly under fatigue. 

The £\3 strain component (fig. 12f) is by far the largest strain component. This strain 
component is due to the eccentricity of the two ends of the tow. The very large magnitude 
suggests that delamination initiation might be dominated by this strain component. 

The strain contours for the case with a waviness ratio of .5 are very similar to those in fig. 12. 
The contours for £ X and £13 are shown in fig. 13 to illustrate the similarity. Interestingly, the 
peak strains tend to be larger for the less wavy weave for the same average axial strain. Table 
2 lists the minimum and maximum strains for the two weaves. However, the less wavy weave 
also has a higher axial stiffness. For a given average strain the average stress is higher for the 
less wavy weave. If both weaves are subjected to the same average stress, the peak strains in 
the more wavy weave are about as large or larger than for the less wavy weave, which is more 
in line with intuition. 

Conclusions 

Techniques were developed and described for performing three-dimensional finite element 
analysis of plain weave composites. The discussion of the analysis emphasized aspects of the 
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analysis which are different from analysis of traditional laminated composites, such as mesh 
generation and representative unit cells. 

The analysis was used to study several different weaves to determine the effects of tow 
waviness on composite moduli, Poisson’s ratios, and internal strain distributions. The average 
normalized composite moduli E Xi G X y, and E z all decreased with increasing waviness. The 
average normalized out-of-plane shear modulus G X z increased with increasing waviness. As 
expected, there are significant strain gradients. The magnitude of these strain concentrations 
suggest that damage initiation will occur at a significantly lower global strain than for a 
traditional cross-ply laminate. However, some of the most severe strain concentrations occur 
where the weave geometry is most difficult to model with confidence. Experimental work is 
needed to characterize the variation of the fiber tow geometry and to document the initiation 
of damage under static and fatigue loads. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November , 1989 
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Appendix A 

Material Property Transformations 

This appendix describes the material property transformations which are required between 
the global and material coordinate systems. Transformations from one coordinate system 
to another are needed for strains, stresses, and the material properties. Although the finite 
element analysis calculates the strains and stresses relative to the global coordinate system, it is 
reasonable to report them relative to the material coordinate system. The material properties 
are known relative to the material coordinate system, but the finite element formulation requires 
the properties relative to the global coordinate system. Hence, there are transformations in 
“both directions”. 

To obtain the strains with respect to the material coordinate axes requires a second order 
tensor transformation. The transformation rule between the global and material coordinate 
systems is given by (ref. 10) 




(Al) 


where the s'- are the strains in the material coordinate system and a,j = cosine of angle 
between the x\ \ and Xj axes which are the material coordinate axes and the global coordinate 
axes, respectively. Fig. Al illustrates the definition of the a t<? for a rotation about the z-axis 
only. 

In expanded form, eqn. (Al) can be written as 


m e 

£ 

£ 

E 

E 

E 

E 

E 

,E 


P 

P 

P 

P 

P 

P 

P 

P 

31 J 


= [T i] 


£22 

£33 

£12 

^21 

£ 23 

£32 

£13 

-£31- 


(A2) 


where 

'aiian 

0l2 a l2 

ai3«13 

011012 


a 21 a 21 

a 2 2«22 

023^23 

a 21 a 22 


a 31 a 31 

a 32 a 32 

O 33 O 33 

O 31 O 32 

Pi] = 

a U a 21 

a 12 a 22 

Ol3«23 

011022 

a 2 ian 

a 22 a 12 

a 23 a 13 

a 21 a l2 


a 21 a 31 

a 22 a 32 

023«33 

021 032 


a 31 a 21 

a 32 fl 22 

033 a 23 

O 31 O 22 


a ll a 31 

a 12 a 32 

013033 

Oil a 32 


- a 31 a ll 

a 32 a 12 

O 33 O 13 

a 31 a i2 


ai 2 an 

012013 

013012 

OllOi3 

013011' 

a 2 2021 

a 22 a 23 

023022 

021023 

023021 

O 32 O 3 I 

032033 

O 33 O 32 

031033 

O 33 O 3 I 

012021 

a l2 a 23 

013022 

011023 

013021 

a 22 a ll 

a 22 a l3 

a 23 a 12 

021013 

023011 

022031 

022033 

023032 

021033 

023031 

032021 

a 32 a 23 

O 33 O 22 

031023 

033021 

a l2 a 31 

O 12 O 33 

O 13 O 32 

011033 

O 13 O 3 I 

O 32 O 11 

032013 

O 33 O 12 

O 31 O 13 

O 33 OH- 


Assuming the strain tensor is symmetric (£{j — £j { ), and performing trivial modification gives 
the transformation in terms of the engineering shear strains, which are 2x the tensor shear 
strains. 


r«n i 
£ ? 2 
£ 33 
2 £j2 
2 £ o3 
L 2£,3 j 


[ T 2} 


' £11 ' 
£22 
£33 
2£j2 
2£23 
- 2 £l3- 


(A3) 
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where 


'an a u 

ai2 a i2 

013013 

2onO]2 

20i2«13 

2auai3 

a 21 a 2 l 

022 a 22 

023023 

2a 2 ia 2 2 

2022023 

2a 2 ia 23 

“31031 

a 32 a 32 

O33O33 

2031 (J32 

2032033 

2031033 

ana 21 

012022 

013023 

an022 + Oi2«21 

012023+ O13O22 

“11023 + 013021 

a 2ia3i 

022 a 32 

023033 

O21O32 + 022031 

a 22 a 33 + O23O32 

“21033 + 023031 

-ana3x 

012032 

013033 

011032 + 012031 

O12O33 + 013032 

“11033 4- 013031. 


Stress is also a second order tensor, so the transformation is the same as for the tensor strains 
(eqns. (Al) and (A2)). Since there is no distinction between engineering and shear stresses, 
the further manipulations resulting in eqn.(A3) are not applicable. 

The material properties are ordinarily defined with respect to the material coordinate axes, 
but the finite element formulation requires that they be defined with respect to the global 
coordinate axes. The material constitutive coefficients comprise a fourth order tensor. Hence, 
the form of the transformation rule is (ref. 10) 


Qjkl ~~ a im a jn a ko a lpC'Tn7iop (A4) 

Note that the a t j in eqn. (A4) are different than those in eqn. (Al), since the transformation 
is in the opposite direction (ie. from the material axes to the global axes). 

Eqn. (A4) is cumbersome. Alternately, we can derive a matrix form of the transformation 
based on the invariance of the strain-energy density. Also, this alternate transformation makes 
use of the already calculated matrix T^. The invariance can be expressed as 


where 


' £11 
€ 22 



2£ 2 3 

L2e 13 J 


and D is a 6 x 6 constitutive matrix. 

Eqns. (3) and (5) can be combined to obtain 

Z r De = £ T T 2 t D%e 

From eqn. (A6) it is obvious that the transformation rule is 


D = T 2 t D'T 2 


(A5) 


(A6) 

(A7) 
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Figure Al. Definition of for rotation about X3 axis. 



Appendix B 
Strain Smoothing 

This appendix describes the procedure used to smooth the strains. The smoothing procedure 
calculates nodal strains based on the strains calculated at the quadrature points. The procedure 
is based on ref. 11. It is assumed that strains are available at 8 quadrature points. The ordering 
of the terms in the matrices which will be described depends on the numbering sequence for 
the corner nodes in a single element and the quadrature points. Fig. B1 shows the numbering 
sequence for the nodes and the quadrature points. For clarity in the sketch, some nodes and 
quadrature points are not shown. However, the numbering patterns for the hidden nodes and 

points are the same as for those shown. 

The smoothing procedure begins with assuming a functional form for the strains within an 
element. The form assumed herein is 


e(xi , X2, 13 ) = S*(xi , X2, 13 )C l sum on i 
where [S*] = [1 x 2 *3 *1*2 *1*3 *2*3 * 1 * 2 * 3 ] 


(Bl) 


Equation (Bl) can be written for the 8 quadrature points in the following form: 
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where A = 1/V3, B = 1/3, and C — A * B 
Equation B2 can be solved for the coefficients C 1 . 


[c‘i = [//rVi 


(B2) 


(B3) 


where II is the coefficient matrix in equation B2. 

With the C* known, the strains at the nodes can be determined using equation (Bl). The 
smoothing procedure can be summarized in terms of a single matrix multiplication. 


[4] = PK1 


(B4) 


where the subscripts n 
[T]is 


and q indicate 


A B 
B C 
C B 
B A 
B C 
C D 
D C 
.C B 


nodal or quadrature values, respectively 
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and the matrix 


(B5) 
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where 


A = 5 + 3*£ 

£> = 5 — 3*1? 

B = -1 - E 

C = —1 + E 

where E = "n/3 

This procedure gives the nodal strains for an individual element. Most nodes in a mesh are 
shared by several elements. A simple averaging procedure was used to smooth discontinuities in 
the strain field. Since strain discontinuities are physically possible at bimaterial interfaces, the 
averaging was performed among elements of the same material type. This averaging procedure 
must be performed before transforming the strains to the material coordinate system. 
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Table 1. Effect of Waviness on Moduli and Poisson’s Ratios 
(Coarse Mesh Results) 


| Waviness Ratio 


.167 

.25 

.5 

Ex 

.92 

.88 

.75 

Ez 

.95 

.93 

.84 

Gxy 

.96 

.94 

.87 

Gyz 

1.10 

1.14 

1.22 

V X y 

.87 

.81 

.60 

Vzy 

1.14 

1.21 

1.45 

U X z 

1.10 

1.14 

1.28 


Table 2. Maximum Normalized Strains for Waviness Ratios of .167 and .5 


Normalized 

Strain* 

Waviness Ratio ~ .167 

Waviness Ratio = .5 

Min. 

Max. 

Min. 

Max. 

£i 

.44 

2.7 

.35 

2.2 

i-i 

-.42 

.29 

-.32 

.27 

£3 

-1.4 

1.6 

-1.4 

1.2 

^12 

-1.0 

1.0 

-.82 

.82 

^23 

-.83 

.83 

-.68 

.68 

^13 

-3.4 

3.4 

-2.9 

2.9 


*Normalized Strain = Strain/Specified Axial Strain 






(a) Symmetric stacking 

Figure 2. Stacking of plain weave mats. 




90° tow 


(b) Unsymmetric stacking 


Figure 2. Concluded. 



(a) Coarse mesh (595 nodes, 96 elements) 


Figure 3. Typical finite element meshes. 
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(a) Moduli 



Figure 8. Effect of waviness on moduli and Poisson’s ratio. 
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90° tow 



Figure 9. Identification of area for which strain contours will be presented. Waviness 




90° tow 



Figure 10. Identification of area for which strain contours will be presented. Waviness ratio = 



Figure 11. Normalized axial strain, E\ /e%. 
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( d ) £ 12 / 4 - 

Figure 12. Continued. 
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(e) £2 sAr- 



(0 fiaAS- 

Figure 12. Concluded. 
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(b) e\$/e%. 

Figure 13. Normalized strain contours for 0° tow. Waviness ratio = .5. 
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